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Abstract 

This paper presents a parametric stiffness analysis of the Orthoghde, a 3-DOF translational ParaUel Kine- 
matic Machine. First, a comphant modehng of the Orthoghde is conducted based on an existing method. Then 
stiffness matrix is symbohcally computed. This allows one to easily study the influence of the geometric design 
parameters on the matrix elements. Critical links are displayed. Cutting forces are then modeled so that static 
displacements of the Orthoglide tool during slot milling are symbolically computed. Influence of the geometric 
design parameters on the static displacements is checked as well. Other machining operations can be modeled. 
This parametric stiffness analysis can be applied to any parallel manipulator for which stiffness is a critical issue. 

Keywords : Parallel Kinematic Machine ; Stiffness ; Parametric analysis ; 

1 Introduction 

Parallel manipulators are claimed to offer good stiffness and accuracy properties, as well as good dynamic 
performances. This makes them attractive for innovative machine-tool structures for high speed machining [1], 
[2], [3]. When a parallel manipulator is intended to become a Parallel Kinematic Machine (PKM), stiffness 
becomes a very important issue in its design [3] , [S] , [5] . In this paper is presented a parametric stiffness analysis 
of the Orthoglide, a 3-axis translational PKM prototype developped at IRCCyN [7J. 

One of the first stiffness analysis methods for parallel mechanisms is based on a kinetostatic modeling [5]. 
It proposes to map the stiffness of parallel mechanisms by taking into account the compliance of the actuated 
joints. This method is not appropriate for PKM whose legs, unlike hexapods, are subject to bending [9 . In 
[TUj . a stiffness estimation of a tripod-based PKM is proposed that solves this problem. However the stiffness 
model presented is not general enough. This gap is filled in [TT] . The method is based on a flexible-link lumped 
parameter model that replaces the compliance of the links by localized virtual compliant joints and rigid links. 
The approach has two differences from that presented in [lOj (i), the way the link compliances are modeled 
(ii), the equations allowing the computation of the stiffness model, that are more general. 

This method is applied to the Orthoglide for a parametric stiffness analysis. The stiffness matrix elements 
are symbolically computed. This allows an easy analysis of the influence of the Orthoglide critical design 
parameters. No numerical computations are conducted until graphical results are generated. First we present 
the Orthoglide, then the compliant model. The results showing the influence of the parameters are presented 
next. 

2 Compliant modeling of the Orthoglide 

2.1 Kinematic architecture of the Orthoglide 

The Orthoglide is a translational 3-axis PKM prototype designed for machining applications. The mobile 
platform is connected to three orthogonal linear drives through three identical RPaR serial chains (where R 
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stands for Revolute and Pa for Parallelogram) (Fig[T|) , and it moves in the Cartesian workspace while maintai- 
ning a fixed orientation. The Orthoglide is optimized for a prescribed workspace with prescribed kinetostatic 
performances. Its kinematic analysis, design and optimization are fully described in [12j . 




Fig. 1 - The Orthoglide kinematic architecture 



2.2 Parameters for compliant modeling 

The parameters used for the compliant modeling of the Orthoglide are presented in Figure [2] and in Table 
[TJ These parameters correspond to a "beam-like" modeling of the Orthoglide legs. The foot has been designed 
to prevent each parallelogram from colliding with the corresponding linear motion guide. 




Fig. 2 - Leg geometrical parameters 



Parameter 


Description 


Lf 
hf, bf 
If, = bf.hj/12 

If, = hf.tyu 

d 

Lb 

Sb 


Foot length, see Fig[5] 

Foot section sides 

Foot section moment of inertia 1 

Foot section moment of inertia 2 

Distance between parallelogram 

bars, see Fig[5] 

Parallelogram bar length, see Fig[5] 
Parellelogram bar cross-section 
area 



Tab. 1 - Geometric parameters 



2.3 Compliant modeling with flexible links 

In the lumped model described in [13j . the leg links are considered as flexible beams and are replaced by 
rigid beams mounted on revolute joints plus torsional springs located at the joints (Fig. [3]). 

From deriving the relationship between the force F and the deformation y{x), the local torsional stiffness 
k can be computed : 

EIy"{x) = F{L-x) 



2 




(a) Flexible beam (b) Virtual rigid beam 



Fig. 3 - General model for flexible link 



EIy{L) = 
^ 6* ~ y{L)/L = FL^/iEI 
k = FL/e 
k = 3EI/L 

If the Orthoglide leg actuator is locked, then one leg can withstand one force F and one torque T (Fig. [J]), 
that are transmitted along the parallelogram bars and the foot. For a compliant modeling that uses virtual 
compliant joints, it is important to understand how external forces are transmitted, and what their effect on 
the leg links is. 




Fig. 4 - Forces transmitted in a leg 



Three virtual compliant joints are modeled along the Orthoglide leg. They are described in Table [H The 
determination of the virtual joint stiffnesses is not detailed here to save space. Since the actuator is assumed to 
be much stiffer than the virtual compliant joints, it is not included in foot stiffnesses. The leg links compliances 
modeled in Table [2] were selected beforehand as the most significant. The joint compliances are not taken into 
account in our model. 

3 Symbolic computing of the Orthoglide stiffness matrix 

The virtual compliant joints stiffnesses depend on the design parameters. In this section, symbolic expres- 
sions of stiffness matrix elements are symbolically computed. This part is based on a stiffness model that was 
fully described in [llj . Therefore, the description of the model will only be summarized here. Fig. [H] represents 
the lumped model of a leg with ffexible links. 

The Jacobian matrix of the ith leg of the Orthoglide is obtained from the Denavit-Hartenberg parameters 
of the ith leg with flexible links. This matrix maps all leg joint rates (including the virtual compliant joints) 
into the generalized velocity of the platform, i.e., 

^idi — t 

where 

= [ ^ii ^ia ^ii ^is ^ie ^is ] 

is the vector containing the 8 actuated, passive and compliant joint rates of leg i, t is the twist of the platform. 
The Pa joint parameterization imposes = — 0^^^.^, which makes and 9i^^^^ dependent. 0ir^ is chosen to 
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Stiffness and description 


Figure 






kact ■ translational stiff- 
ness of the prismatic ac- 
tuator 






r 


2EIf^/Lf : foot bending 
due to torque T 








3EIf-^/Lf : foot bending 
due to force F 




ESBd^/2LB ■■ differen- 
tial tension of parallelo- 
gram bars due to torque 
T 





Tab. 2 - Virtual compliant joints modeling 




Fig. 5 - Flexible leg 



model the circular translational motion, and finally Jj is written as 

J ^12 *^'i3 



S'ig X fig ^isbis X Tjgj^^ e^g X Tig e,;^ X Yi^ Big X r^g 

in which Si- is the unit vector along joint j of leg i and Vi. is the vector connecting joint j of leg i to the platform 
reference point. Therefore the Jacobian matrix of the Orthoglide can be written as : 



Jl 






















Js 



One then has : 

je = Rt (1) 
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with R : 



h 
h 
h 



and t 



V 



9 being the vector of the 24 joint rates, that is 9 — [ Oj^ ^3" ] • Unactuated joints are then ehminated 

by writing the geometric conditions that constrain the two independent closed-loop kinematic chains of the 
Orthoglide kinematic structure : 



Jl^l — J202 

Ji^i = ^383 

From ^ and ^ one can obtain (see [TT] for details) : 



(2) 
(3) 



A9' = B9" 

where 9' is the vector of joint rates without passive joints and 9" is the vector of joint rates with only passive 
joints. Hence : 

9" = B-^A9' 

Then a matrix V is obtained (see [TT] for details) such that : 



From (H]) and (01) one can obtain 



JV0' = Rt 



(4) 



(5) 



As matrix R represents a system of 18 compatible linear equations in 6 unknowns, one can use the least-square 
solution to obtain an exact solution from ^ : 



t = (R^R)-iR'^JV^' 



Now let J' be represented as : 



J' = (R^ R)^R^ JV 



Then one has : 



t = 3'9' (6) 

According to the principle of virtual work, one has : 

T^9' = w^t (7) 

where r is the vector of forces and torques applied at each actuated or virtual compliant joint and w is the 
external wrench applied at the end effector, point P. Gravitational forces are neglected. By substuting ([6]) in 
([7]), one can obtain : 

r = 3'^^^ (8) 

The forces and displacements of each actuated or virtual compliant joint can be related by Hooke's law, that 
is for the whole structure one has : 



T = KjA9' (9) 

with : 



A 











A 











A 
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and : 



kact 

















2EIf 







ESBd^ 



(10) 



/S.6' only includes the actuated and virtual compliant joints, that is by equating ([S]) with @ 



IT, 



hence 



K 7A6I' = J 



Pre-multiplying both sides by J' one obtains : 

J'A0' = J'KjiJ'^w 



(11) 



Substituting ^ into pTjl . one obtains 



t = J'KT^J'^w 



Finally the compliance matrix k, is obtained as follows : 

n = J'Kj-ij'^ 

In the Orthoglide case we obtain a simplified compliance matrix : 



/ 













Kl5 


^ 







Kll 





'«24 





'«26 















'tss 















K44 









'tis 





K35 





'tss 





V 


'tis 












t66 / 


is 


















K = 


J'Kj 


-lj,T 


-1 





(12) 



Remarks : 

From the expression of matrix A (see eq. [TU)) we deduce that all elements of stiffness matrix K can be 
factored either in factors of E or in factors of kact- As kact is assumed infinite compared to other stiffnesses, 
the elements of K which are factors of kact must be eliminated. Also obviously when E increases, K elements 
increase too, which is in accordance with intuition. Young's modulus E should then be eliminated from K 
elements expressions because its influence needs no further study. 

In eq. (|12p one should notice that kh = K22 = 'tsa, and K15 = Kig. Also, we noticed that when X=Y=Z we 
have K44 — K55 — Kge. 



4 Parametric stiffness analysis 

In this part, we study the in uence of the geometric parameters on the elements of the Orthoglide stiffness 
matrix. For a realistic inspection of the symbolic expressions of these elements, it is necessary to replace some 
parameters by numerical values, once the symbolic expressions have been obtained. Numerical values of the 
design parameters come from the Orthoglide prototype [Hj. 
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4.1 Qualitative analysis at the isotropic configuration 



There is an isotropic configuration in the Orthoghde workspace [7] . This configuration gives a good insight 
of the Orthoghde general performances. It happens when the tool point P is at the intersection of the three 
actuated joint axes (see Fig. 1). In this configuration, the stiffness matrix K is diagonal and the symbolic 
expressions of its elements are simplified : 



K = E/2 



( Ka 1 

Ka 
1 
1 
1 
1 



V 1 



E = 2Ka being the torsional stiffness with 



Ka 



cPSe 



1 



1 



1 



LB 'iLfcos^X 



hfb) 



and E — 2Kb being the translational stiffness with 



1 \ 



1111 

Ka I I I 

I Kb I I 
I I Kb I 
1 1 I Kb J 



(13) 



Kh 



0.5hjbf 



E 



(14) 



From the inspection of Ka we note that increasing the cross-section areas SB and bfhf increases the 
torsional stiffness, which is in accordance with intuition. Increasing the distance d between the paralle- 
logram bars increases the torsional stiffness. Also, it is no surprise that increasing Lf or Lb decreases 
the torsional stiffness. It is no surprise that Ka and Kb increase when E increases, and that Kb increases 
when kact increases, while Ka does not depend on fcact, the prismatic actuated joint stiffness. Ka and 
Kb do not depend on e at the isotropic configuration. 

From the inspection of Kb we note that increasing the cross-section area bfhf increases the translational 
stiffness, as increasing the foot length Lf does. We also note that if A (angle between the foot and the 
actuated joint axis) decreases toward zero, then the translational stiffness tends toward infinite. This is 
in accordance with intuition since at the isotropic configuration, if A = then the virtual compliant joint 
"foot bending due to force F" (see Tab. 2) is not requested because the applied force axis crosses the 
virtual compliant joint axis. 

Furthermore, if A = 0, the cantilever disappears. As this cantilever is aimed at preventing each paral- 
lelogram from colliding with the corresponding linear motion guide, from a designer point of view the 
conclusion could be : "A should be lowered to increase the translational stiffness, but there should remain 
enough cantilever for the workspace volume to be large enough" . 



4.2 Quantitative analysis at the isotropic configuration 

The qualitative analysis conducted above provides an interesting insight of the influence of the geometrical 
parameters on rotational and translational stiffnesses. Quantitative information is now obtained by studying 
the influence on Ka and Kb of a +/ — 50% variation of the numerical value of each design parameter. The 
results are shown on FigIH]and on FiglT) From Fig. 6 analysis, we can go further than the qualitative analysis 
and see that d, Sb and Lb have very little influence on the rotational stiffness while Lf , hf bf and A have 
much more influence. From Fig. 7 we note, as deduced from the qualitative analysis, that decreasing Lf or X 
increases the translational stiffness, and that increasing hf can improve the translational stiffness. Increasing hf 
does not reduce the workspace volume like reducing Lf does. However increasing hf increases the foot weight, 
which could be prejudicial to the dynamic performances. The proof of this forecast needs further investigations. 
Quantitative analysis provides another point of view and helps in a better understanding of the influence of 
the geometric design parameters on the machine stiffness. 



5 Conclusions 

A parametric stiffness analysis of a 3-axis PKM prototype, the Orthoghde, was conducted. First, a compliant 
model of the Orthoglide was obtained, then an existing stiffness analysis method for parallel manipulators 
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Parameter variation (%) Parameter variation (%) 

(a) Influence of L/ , 6/ , hf and A (b) Influence of rf, Lb and Sb 



Fig. 6 - Quantitative influence of the geometric parameters on the rotational stiffness E = 2Ka 




Parameter variation (%) 

Fig. 7 - Quantitative influence of the geometric parameters hf , bf , Lf and d on the translational stiffness 
E = 2Kb 

was applied to the Orthoglide. The stiffness matrix elements were computed symbolically. In the isotropic 
configuration, the influence of the geometric design parameters on the rotational and translational stiffnesses 
was studied through qualitative and quantitative analyses. These provided relevant information for stiffness- 
oriented design or optimization of the Orthoglide. 

This example of a parametric stiffness analysis shows that simple symbolic expressions carefully built and 
interpreted provide broad information on parallel manipulators stiffness features. Comparison of the stiffness 
model with a finite elements model and with stiffness tests on the prototype are currently being conducted. 
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